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Abstract. We extend our previous work on Poisson-like formulas for 
subresultants in roots to the case of polynomials with multiple roots 
. in both the univariate and multivariate case, and also explore some 

closed formulas in roots for univariate polynomials in this multiple roots 
, setting. 

(N 

> 
O 

^ I 1. Introduction 

•/^ . In [DKS2006] we presented Poisson-like formulas for multivariate subre- 

sultants in terms of the roots of the system given by all but one of the input 
■ polynomials, provided that all the roots were simple, i.e. that the ideal gen- 

i erated by these polynomials is zero-dimensional and radical. Multivariate 

resultants were mainly introduced by Macaulay in |Macl902] . after ear- 
lier work by Euler, Sylvester and Cayley, while multivariate subresultants 
were first defined by Gonzalez- Vega in |GLV1990l IGLV199T] . generalizing 
Habicht's method |Habl948] . The notion of subresultants that we use in 
this text was introduced by Chardin in [Chal995] . 
<^ : Later on, in |DHKS2007l IDHKS2009j . we focused on the classical uni- 

1^ ■ variate case and reworked the relation between subresultants and double 

Sylvester sums, always in the simple roots case (where double sums are 
^2 I actually well-defined). This is also the subject of the more recent articles 

IRS20111 IKS2012j . As one of the referees of the MEGA'2007 conference 
, pointed out to us, working out these results for the case of polynomials with 

I multiple roots would also be interesting. 

This paper is a first attempt in that direction. We succeed in describing 
Poisson like formulas for univariate and multivariate subresultants in the 
presence of multiple roots, as well as to obtain formulas in roots in the 
univariate setting for subresultants of degree 1 and of degree immediately 
5^ . below the minimum of the degrees of the input polynomials: the two non- 

trivial extremal cases in the sequence of subresultants. We cannot generalize 
these formulas for other intermediate degrees, and it is still not clear for us 
which is the correct way of generalizing Sylvester double sums in the multiple 
roots case. 

The paper is organized as follows: In Section [2] we recall the definitions 
of the classical univariate subresultants and Sylvester double sums, and of 
the generalized Wronskian and Vandermonde matrices. We then show how 
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the Poisson formulas obtained in [Honl999] for the subresultants in the 
case of simple roots extend to the multiple roots setting by means of these 
generalized matrices. We also obtain formulas in roots for subresultants in 
the two extremal non-trivial cases mentioned above. In Section [3] we present 
Poisson-like formulas for multivariate subresultants in the case of multiple 
roots, generalizing our previous results described in |DKS20()6] . 

Acknowledgements: We wish to thank the referee for her/his careful read- 
ing and comments. A preliminary version of these results was presented at 
the MEGA 2009 Conference in Barcelona. Part of this work was done at the 
Fields Institute in Toronto while the authors were participating in the Fall 
2009 Thematic Program on Foundations of Computational Mathematics. 



2. Univariate Case: Subresultants in multiple roots 

2.1. Notation. We first establish a notation that will make the presentation 
of the problem and the state of the art simpler. 

Set d,e £ N and let A := (ai,...,ad) and B := . . . , /3e) be two 
(ordered) sets of d and e different indeterminates respectively. 

For m,n G N, set {di,...,dm) G and (ei,...,e„) G N" such that 
di + ■ ■ ■ + dm = d and ei + • • • + e„ = e, and let 

A:= {{ai,di);...;{am,dm)) and 5 := ei); . . . ; (/3„, e„)) 

(these will be regarded as "limit sets" of A and B when roots are packed 
following the corresponding multiplicity patterns). 

We associate to A and B the monic polynomials / and g of degrees d and 
e respectively, and the set R{A,B), where 

d e 

/(x) := JJ(x - Oi) and g{x) := - 
i=i j=i 

R{A,B)= n («*-/3j)'= n 

l<i<d,l<j<e l<i<d 

with natural limits when the roots are packed 

m n 

7{x):=l[{x-aif^ and g{x) := l{{x - (3^)'^ , 
1=1 j=i 

R(A,B)= n («.-/3,)*^^= n ^M"'- 

l<«<m,l<j<n l<i<m 

2.2. Subresultants and Sylvester double sums. We recall that for < 
t<d<e or 0<t<d = e, the t-th. subresultant of the polynomials 
/ = a^x'^ + • • • + ao and g = b^x^ + • • • + 6o) introduced by J.J. Sylvester in 
|Sylvl853] , is defined as 
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d+e-2t 



Siest{f,g) := det 



ad ■ ■ 




at+i-{e-t~ 


-1) 






ad ■ ■ 


at+1 




x^f{x) 


be •• 




• ^t+l-(d-t- 


1) 


x'^~^~^g{x) 




be ■■ 


bt+i 




x^g{x) 



e-t 



d-t 



with = = for ^ < 0. When t = we have Srest{f,g) = Res{f,g). 

In the same article Sylvester also introduced for < p < d,0 < q < e the 
following double-sum expression in A and B, 

Svlv^'^f A B- X) - V R(x A') R(x B') ^i^' ^B') R{A\A' ,B\B') 
Sylv {A,B,x).- 2^ R{A',A\A')R{B',B\B'y 



A'CA,B'CB 
\A'\=p,\B'\=q 



where by convention R{A', B') = 1 if ^' = or B' = 0. For instance 
(1) Sylv'''\A,B;x) = RiA,B)= J] (a^ -/?,) = Res(/, 5). 

l<i<d,l<j<e 

We note that Sylv^''^(A, B; x) only makes sense when Ui 7^ aj and /3j 7^ /3j for 
i ^ j, since otherwise some denominators in Sylv^'''(A, B; x) would vanish. 

The following relation between these double sums and the subresultants 
(for monic polynomials with simple roots / and g) was described by Sylvester: 
for any choice oi < p < d and < g < e such that t := p + q satisfies 
t<d<eoTt = d<e, one has 



(2) 



Srest(/,5) = (-l) 



pid-t) 



-1 



SylvP'i{A,B;x). 



This gives an expression for the subresultant in terms of the differences of 
the roots — generalizing the well-known formula ([1]) — in case f and g have 
only simple roots. However, when the roots are packed, i.e. when we deal 
with A and B, the expression for the resultant is stable, i.e. 



Res(/,5r) 



n 

l<i<m, l<j<n 



/3. 



while not only there is no simple expression of what Srest(/,^) is in terms 
of differences of roots but moreover there is no simple definition of what 
Sylv^'^(^, B; x) should be in order to preserve Identity ([2]). Of course, since 
Sresi(/,g) is defined anyway, Sylv^''^{A, B; x) could be defined as the result 



SylvP'<?(A 



,B;x) := 



Srest (/,£/) 



but this is not quite satisfactory because on one hand this does not clarify 
how Sresf behaves in terms of the roots when these are packed, and on the 
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other hand, Sylv^''^{A, B; x) is defined for every < p < d and < g < e 
while Srest is only defined for t : = p + q < min{d, e}. 

In what follows we express some particular cases of the subresultant of 
two univariate polynomials in terms of the roots of the polynomials, when 
these polynomials have multiple roots. These are partial answers to the 
questions raised above, since we were not able to give a right expression 
for what the Sylvester double sums should be, even in the particular cases 
we could consider. Nevertheless the results we obtained give a hint of how 
complex it can be to give complete general answers, at least in terms of 
double or multiple sums, see Theorem 12.71 below. 

2.3. Generalized Vandermonde and Wronskian matrices. We need 
to recall some facts on generalized Vandermonde and Wronskian matrices. 

Notation 2.1. Setu G N. T/ie generalized Vandermonde or confluent (non- 
necessarily square) uxd matrix Vu{A) associated to A = ((ai,(ii); . . . ; {am, dm)), 
|Kall984j ■ is 



Vu{A) = Vu{{ai,di); . . . ; {am, dm)) ■= 



Vu{ai,di) 



Vu{a 



m, dm) 



where 





1 















1 








Vu{ai,di) := 


a2 


2ai 


1 









{u - l)ar' 







with the convention that when k < j, 



k-j 



0. 



When di = 1 for all i, this gives the usual Vandermonde matrix Vu{A). 
When u = d, we omit the sub-index u and write V{A) and V{A). 

For example 



y((a,3);(/3,2)) 



1 








1 





a 


1 





/? 


1 


a2 


2a 


1 




2(3 


a^ 


3a2 


3a 


/33 


3/32 




4a^ 






4/33 



and 



y3((a,3);(/3,2)) 



1 








1 





a 


1 





/? 


1 


a' 


2a 


1 


/32 


2/3 



The determinant of a square confluent matrix is non-zero, and satisfies, 
[Ait 1 939 j . 
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det{V{A))= J] {aj-a.f^''^. 

In the same way that the usual Vandermonde matrix V{A) is related to 
the Lagrange Interpolation Problem on A, the generalized Vandermonde 
matrix V{A) is associated with the Hermite Interpolation Problem on A 
|Kall984] : Given {j/jj. ,1 < i < m, < jj < dj}, there exists a unique 
polynomial p of degree deg(p) < d which satisfies the following conditions: 

... , p(^i-i)(ai) = (di-l)!j/i,rf,_i, 



p{ai) = 0! yifi, p'{ai) = 1! 
p{a„i) = 0! y,„,o, p'{ajn) = l! J/m,!, 



This Hermite polynomial p = oq + aix + • • • + a(i_ix°' ^ is given by the only 
solution of 

(ao ai . . . ad-i) ■ V{A) = {yi^ yi^i . . . ym,dm-i) 

(here the right vector is indexed by the pairs (i, ji) for 1 < i < m, < ji < d,i) 
and satisfies 



(3) det [V{A)) p{x) = - det 



d 


1 




1 




X 


V{A) 






yi,o 2/1,1 • • • ym,dm-i 






The polynomial p can also be viewed in a more suitable basis, where the 
corresponding "Vandermonde" matrix has more structure. We introduce 
the d polynomials in this basis. 

Notation 2.2. For 1 < i < m we set 



and, for < ki < di, 



f 



{x — ai)'^^ 



Then, in this basis, the polynomial p = Yli ki '^iMfiM is given by the only 
solution of 



(ai,o ai,i . . . am,d™-i) • V'{A) = {yi^o . . . ?/m,d™-i) 



where 



di 



V'{A) :- 



V'{ai,di) 






















V'{am,dm) 
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with 



V'{ai,di) :-- 






7'i{ai) ■ ■ 

7i{m) ■ ■ 




(di-iy. 
















7i(tti) 



and satisfies 


d 


1 






/l,0 








det {V'(A)) p{x) = - det 


V'(A) 








/ m,dm-l 




2/1,0 yi,l ■■■ ym,dm-i 






We note that 

det(y'(:4))= n 7iK)* 



l<i<m 



m(m—l) 



7 , \ 2 m(m— 1) r, 

la, - ) = (-1)^^ det(V(A))2. 



( n (< 

l<i<j<m 

In particular p = j. UijiPiji where for 1 < i < m, < ji < di, the 6as«c 
Hermite polynomials pij. are the unique polynomials of degree deg(pjj.) < d 
determined by the conditions for 1 < i < m, < qi < dg, 



(4) 



Pij^ {(^d = ji- for e = i and qe=ji, 
p'fj^{ai) = otherwise. 



When A = A = (ai, . . . , a^), then, denoting /j := /j, we have 
X — aj^ 1 



Pi,o 



/i for 1 < ^ < 



while for A = {a, d), 

Pi,j = {x- ay = Jij for Q<j<d. 
The following proposition generalizes these two extremal formulas. 
Proposition 2.3. Fix 1 <i <m and < j < di. Then 



di-l-j 



Pi,j 



E 



fdi-l+ke\ 



\ki+...+ki+...+k^=k 



f, 



where ki + ■ ■ ■ + ki + ■ ■ ■ + km denotes the sum without ki. (When m = 1, 
the right expression under brackets is understood to equal 1 for k = and 
otherwise.) 
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Proof. Applying for instance |Spl960 Th. 1], we first remark that 



di-l-j 

E 

k=0 



1 / 1 



(k) 



Then we plug into the expression the following, given by Leibnitz rule for 
the derivative of a product: 



/ 1 \ (fc) 



kl+... + ki + ... + kr = k 



^de-l+ki^^ 



□ 



The basic Hermite polynomials enable us to compute the inverse of the 
confluent matrix V{A): 



V{A) 



V* 



where V* : = 



coefficients of pi^i 
coefficients of pi^di 



di 1 



l<i<r 



(here the coefficients of j. (x) are written in the monomial basis 1, x, . . . x*^ ^). 

Now we set the notation for a slight modification of a case of generalized 
Wronskian matrices. 

Notation 2.4. Set n € N. Given a polynomial h{z), the generalized Wron- 
skian (non-necessarily square) u x d matrix Wh,u{A) associated to A = 
((ai,di); . . . ; (am, dm)) is 



Wh^uiA) = Wh,u{{ai,di); . . . ; (am, dm)) := 




where 



di 



Wh,u{ai,di) := 


h{ai) h'{ai) 
{zh){ai) {zh)'{ai) 






(z"-i/i)(a,) {z^'^hyiai) .. 


(d,~l)\ 



When u = d, we omit the sub-index u and write Wh{A). 
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For example for h{z) = x — z and A = (a, 3), 



3 


X — a 




-1 





ax — o? 


X - 


- 2a 


-1 


1 3 
ax — a 


2ax 


-3a2 


X — 3a 



The determinant of a square Wronskian matrix is easily obtainable perform- 
ing row operations in the case of one block, and by induction in the size of 
the matrix in general: 



det {WhiA)) 



n ("^- 

l<i<j<r?i 



)''^^^)h{ai)^'---h{amr 



2.4. Subresultants in multiple roots. In this section, we describe ex- 
plicit formulas we can get for the non-trivial extremal cases of subresul- 
tants in terms of both sets of roots of / = (x — aiY^ ■ ■ ■ {x — am)'^"' and 
g = (x - (3i)^^ ■ ■ ■ {x - pnY- with d = Y7=i di and e = E"=i e^- More 
precisely, we present formulas for Svest(f,g) for the cases t = d — 1 < e 
(Proposition 12.61 below) and t = 1 < d < e (Theorem 12. 7p . We will de- 
rive them from Theorem 12.51 below, a generalization of |Honl999[ Th. 3.1] 
and |DHKS2007l Lem. 2] which includes the multiple roots case (and is also 
strongly related to a multiple roots case version of [ DHKS20091 Th. 1]). The 
main drawback of this approach to obtain formulas for all cases of t is the 
fact that submatrices of generalized Vandermonde matrices are not always 
generalized Vandermonde matrices, so in general their determinants cannot 
be expressed as products of differences. This is why the search for nice 
formulas in double sums in the case of multiple roots is more challenging. 

Theorem 2.5. SetO<t<d<eorO<t<d = e. Then 



Srest(/,5) = (-l)'*"* det {V{A)) ' det 



(-1)'^ det {V{A)) ^ det {V{B)) ^ det 



d 


1 






1 




Vt+i(A) 




t+i 








Wg,d-t{A) 





d-t 



d 


e 


1 






1 


Vt+i(A) 





X* 






Vd+e-t(A) 


Vd+e-t(B) 






t+l 



d+e-t 



where c := max{e (mod 2),d — t (mod 2)}. 
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Proof. The proof is quite similar to the proofs of Lemmas 2 and 3 in [DHKS2007] , 
replacing the usual Vandermonde and Wronskian matrices by their general- 
ized counterparts. We will thus omit the intermediate computations. 
Let / = Yli=(}'^i^^j where Od = 1, and g = Yl'j=o^j^i^ where be = 1- We 
introduce the following matrices of [DHKS2007J : 



Mj:= 



and 



St := 



d+e-t 


aa ... 


ad 






ao . 


■ ■ ad 


d+e-t 








t 




Mj 


e-t 


where 




d-t 





d+e-t 


bo ■ 


.. be 






bo ■■ 


be 



d+e-t 



X -1 




. 


X 


-10.. 


. 



We have ([DHKS20021 Lem. 1]): 

Sresi(7,5) = (-1)(^"*)("-*) det(5t). 
Also, exactly as in the proof of |DHKS2007i Lem. 2], 



d-^e—t d e—t 



t 









d 


W^-,.t{A) 


* 


e-t 


Mj 


Vd+e-t{A) 











d-t 






Ide-t 


e-t 


Wg,d-t{A) 





This implies first the generalization of |DHKS2007[ Lem. 2] to the multiple 
roots case: 



det{V{A))STest{f,g) = det 



W^.,,tiA) 

Wg,d-t{A) 



t 

d-t 



(-l)''-'det 



d 


1 






1 




Vt+i(A) 




t+i 








Wg,d-t{A) 





d-t 



where the second equality is a consequence of obvious row and column 
operations. Next, we get as in the proof of [DHKS20071 Lem. 3], 





d 


e 


1 










1 




det (ViA)) det (ViB)) Srest(/,5) = (-1)^ det 


Vt+i(A) 







t+1 




Vd+e-tiA) 


Vd+e-t{B) 





d+e-t 



□ 
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We note that starting from the first equahty above and applying similar 
arguments, we also get very simply 
(5) 

d e 

t 



det {V{A)) det {V{B)) Srest{f,g) = det 










Vd+e-t{B) 



d+e-t 



As mentioned above, when i = the formula in roots for SresQ{f,g) 
specializes well when considering Sreso(/,ff)- When t = d < e, the for- 
mula Sresrf(/, (7) = ni<i<(i(^ ~ Ofj) also specializes well as Sresrf(/,5) = 
ni<i<m(^ — ai)'^\ Our purpose now is to understand formulas in roots for 
the following extremal subresultants, i.e for Sresi and Sresrf_i, in case of 
multiple roots. 

• The case t = d — 1 < e: When / has simple roots, it is known (or can 
easily be derived for instance from Sylvester's Identity ([2]) for p = d — 1 and 
q = 0) that 

d 

n- — —) = y"9(ai)Pi, 



Sicesd-i{f,g) 



d 

E 



where pi is the basic Lagrange interpolation polynomial of degree strictly 
bounded by d such that Pi{ai) = 1 and Pi{aj) = for j ^ i. In other words, 
Sresrf_i(/, is the Lagrange interpolation polynomial of degree strictly 
bounded by d which coincides with 'g in the d values ai, . . . , ad- This for- 
mula does not apply when / has multiple roots, but we can show that we 
get the natural generalization of this fact, that is, that Svesd-i{f ,g) is the 
Hermite interpolation polynomial of degree strictly bounded by d which co- 
incides with g and its derivatives up to the corresponding orders in the m 
values am- 



Proposition 2.6. 



Sresd„i(/,5r) = 2^ 2^ — Pi,],, 



i=l ji=0 

where pij. is the basic Hermite interpolation polynomial defined by Condi- 
tion ^ or Proposition [Ql for A. 

Proof. In this case, applying the first statement of Theorem 12.51 we get 

d 1 



Sresd-i(/,5) = -det(y(A)) ' det 



Vd{A) 


1 

^d— 1 


Wg,l{A) 






where when following the subindex notation of Formula ([3]), we note that 
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The conclusion follows by Formula 
For example, when A = (a, d), we get 



11 

n 



Sres,_i((x - ay,9) = ^ ^^-^ {x - aY 



j=0 



the Taylor expansion of g up to order d — 1. 

• The case t = 1 < d < e: We keep Notation 12.21 When / has simple roots, 
it is known (or can easily be derived for instance from Sylvester's Identity ([5]) 
for p = 1 and q = 0) that 

(6) Sresi(/,5) = (-ir^X:(n^ 

1=1 

The general situation is a bit less obvious, but in any case we can get an 
expression of Sresi {f,g) by using the coefficients of the Hermite interpolation 
polynomial, in this case of the whole data 

AUB := ((ai,di); . . . ; (arri,dm); (/?i,ei); . . . ; (/3n,e„)). 

We note that 

det {V(AUB)) = det {V(A)) det {V(B))R(B,A) 
which holds even when = (3j for some 
Theorem 2.7. 



i=l 



E 



n 



fci +... + fe; +... l<J<m ^ - l<e<n 



■ + fc„+„ = di - 1 j i 



min{l,di-l} ^ 
■ ■ + 5 

11 + 71 ■ 

Proof. Setting t = 1 in Expression ([5]) we get 



/dj-l+kj\ 

TT ^ ''^ ' TT 



H h Si H 1 < j < m 

■ ■ + fe,„ + „ = di -2 j 7^ i 



det det (^(5)) Sresi(/,5) = det 



VFx-.,i(^) 





Vd+e-l{A) 


V^d+e-l(5) 



1 

d+e- 



where 

Wx-z,i = (x - ai, -1,0, . . . ,0, . . . ,x - am, -1,0. . . ,0). 

^ V ' ^ V ' 

di dm 

We expand the determinant w.r.t. the first row, and observe that when we 
delete the first row and column j, the matrix that survives coincides with 
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V{A U -B)(rf_|_e,j)i the submatrix of V{A U B) obtained by deleting the last 
row and column j. Therefore, 








Vd+e-l{A) 





det 



= (det {V{A U B)\^d+e,m))i^ - + det {V{A U i3)|(rf+e,0'O-)) 

i=i 

where (/>(i) equals the number of the column corresponding to (1, a^, ... , af~^'^~^) 
in V(A\J^), and = 0(z) + 1 if dj > 1 and otherwise. 
Now, from 

det {V(AuB)\(^a+eMj))) = det {V(AuB)) ViAuB)-^^^^^^^, 

det {V(AuB)\(^a+e,<t>'U))) = {-1)"+'-^'^'^ det {V(AUB)) V(AUB)~}^.^^^^^ 
(by the cofactor expression for the inverse) and 

det {V(AUB)) = (-l)'^'^det {V(A)) det {V(B))R(A,B) 
we first get, since R{A,B) = ni<i<m ^('^O'^S that 

(rn m 
_ _ i—l 2—1 

We set h := fg, and for z = 1, . . . , m, /ij := h/{x — aiY^ In |Csl975t Id. 9], 
it is shown that 

(4-1) 

(Oi), 



+e 



and when > 1, 

f(au:b)-i 



(di - 1)! \h 



1 



1 



1 



K-2) 



-/-'(O.'^+e ((i- _ 2)! 

Therefore, we obtain the statement by applying Leibnitz rule 

(k) fdj-l+kj\ (ei~l + k. 

E n TTTT^w n 



□ 

Note that in the case that / has simple roots we immediately recover 
Identity ^ while when f = {x — a)'^ for d > 2, we recover Proposition 3.2 
of [DKS2009j: 



E n 



VJ-J- (a- Be) 

fci + ---+fc„=ci-2 £=1 ^ 



^ (x - a) + 
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3. Multivariate Case: Poisson-like formulas for Subresultants 

We turn to the multivariate case, considering the definition of subresul- 
tants introduced in |Chal995] . Our goal is to generalize Theorem 3.2 in 
|DKS2006] -that we recall below- to the case when the considered polyno- 
mials have multiple roots. We first fix the notation, referring the reader to 
|DKS2nn6j for more details. 

3.1. Notation. Fix n G N and set A G N for 1 < -i < n -h 1. Let 

fi • — ^ ^ Oii,a'^ ^ K[xi , . . . , Xfi\ , 
\a\<D, 

be polynomials of degree Di in n variables, where cx = (ai,...,a„) € 
(Z>o)", a;" := x"^ • • • x"", |q;| = ai -|- • • • -|- a„, and K is a field of charac- 
teristic zero, that we assume without loss of generality to be algebraically 
closed. 

Fix t €N. Let k := 'HDi...D„+i{t) be the Hilbert function at t of a regular 
sequence of n + 1 homogeneous polynomials in n + 1 variables of degrees 
Di,. . . ,Dn+i, i-e. 

k = #{ic" : \a\ < t, Ui < Di, 1 < i < n, and t — |q;| < Dn+i}- 

We set 

S := {a;i'i,...,a;^fc} C K[x]t 
a set of k monomials of degree bounded by t, and 

A5(/l, • • • , /n+l) := '-^^{fl, fn+l), 

for the order t subresultant of f^, . . . , f^^i with respect to the family S'^ := 

{x'^'^xI^_IJ^ \ • • • ) 33'"''= x^"^7'° ' } defined in |Chal995| . Here, f^ denotes the ho- 
mogenization of fi by the variable Xn+i- 

We recall that the subresultant A5 is a polynomial in the coefficients of 
the fi of degree 'HDi...Di_iA+i...£'n+i - A) for 1 < « < n + 1, having the 
following property: A5 = if and only if ItUS^ does not generate the space 
of all forms of degree t in k[xi, . . . , where It denotes the degree t part 

of the ideal generated by the //*'s. 

By |Chal994] we know that 
(7) det{Ms) = £{t)As, 

where M5 denotes the Macaulay-Chardin matrix obtained from 



(8) 



by deleting the columns indexed by the monomials in S, and £{t) is the 
extraneous factor defined as the determinant of a specific square submatrix 
of dSD (see |Chal995[ [Chal994l[DK52006] ). 

We set p:= iDi-l) + --- + {Dn- 1) and for j > 0, tj := nD^...Djj), the 
Hilbert function at j of a regular sequence of n homogeneous polynomials 
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in n variables of degrees Di, . . . , Dn- We define also 



(9) T, 



any set of Tj monomials of degree j if j > max{0, t — Dn+i + f } 
{a;" : \a.\ = j, < A for f < i < n} if < j < t - Ai+i + 1, 



and denote with D := Di ■ ■ ■ Dn the Bezout number, the number of common 
solutions in of n generic polynomials. 

Set r := Uj>oTj and T* := U^^^,^^?^-. Note that \T\ = D. We enumerate 
the elements of T as follows: T = {ic"i, . . . , a;"^}, and assume that for 
s := |T*| we have T* = {a;"i, . . . , tc""} . Also set 
(10) 

n := . . . ,ic^'} = : \a\ < t,ai < Di,! <i <n, t- \cx\ > -D„+i}. 



Finally, for 1 < i < n, let fi be the homogeneous component of degree D 
of fi, and Aj- := A^r {fi, . . . , fn) be the order j subresultant of /i, . . . , /r 
with respect to Tj. 



3.2. Poisson-like formula for subresultants. From now on we assume 
that /i , • • • , /n are generic in the sense they have no roots at infinity (which 
implies by Bezout theorem that the quotient algebra A := K[x]/{fi, . . . , fn) 
is a finitely dimensional i^T- vector space of dimension D, which equals the 
number of common roots in K" of these polynomials, counted with mul- 
tiplicity, see e.g. |CL01998l Ch. 3, Th. 5.5]), and that T is a basis of 
A. 

In |DKS2006] we treated the case of general polynomials with indetermi- 
nate coefficients, which specializes well under our assumptions to the case 
when the common roots ^i, . . . , of /i, . . . , /„ in are all simple. Set 
Z := {^i, . . . , We introduced the Vandermonde matrix 



(11) 



Vr{Z) :-- 









tOLD 



G K 



DxD 



whose determinant is non zero, since T is assumed to be a basis of A, and 
we defined 



(12) 



Os{Z) 





































/ii+i 


(Ci) • 






fn+l 


(ll) • 


■ fn+liio) 



£ K 



DxD 
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Theorem 3.1. |DKS2006l Th. 3.2] For any t G Z>o and for any S 

{x'^^, . . . , x'^''} C -fC[a;]t of cardinality k = 'HDi...D„+i(t)> we have 



A5(/i,...,/n+i) = ± I n ^r, 

J=t-Dn + 1 + 1 



det {OsjZ)) 
det {Vr{Z)) • 



In order to generalize this result to systems with multiple roots, and 
obtain an expression for the subresultant in terms of the roots of the first 
n polynomials /i, • • • we need to introduce notions of the multiplicity 
structure of the roots that are sufficient to define (/i, • • • , /n)- To be more 
precise, in the case of multiple roots, the set of evaluation maps {ev^ : A 
K I ^ common root of /i, ...,/„} is not anymore a basis of A* , the dual of 
the quotient ring A as a ET-vector space, though still linearly independent. 
Hence other forms must be considered in order to describe A* and to get a 
non-singular matrix generalizing Vf{Z). 

All along this section we will use the language of dual algebras to gener- 
alize Theorem 13.11 for the multiple roots case (see for instance in [KK1 9871 
IBCRS1996j and the references therein). In Theorem 13.41 below we show that 
any basis of the dual A* gives rise to generalizations of Theorem 13.11 as long 
as we assume that T is a basis of A. This is the most general setting where 
a generalization of Theorem 13.11 will hold. However, this version of the The- 
orem, using general elements of the dual, does not give a formula for the 
subresultant in terms of the roots. 

In order to obtain these expressions, we need to consider a specific basis 
of A* which contains the evaluation maps described above. It turns out 
that one can define a basis for A* in terms of linear combinations of higher 
order derivative operators evaluated at roots of /i, . . . , /„. This is the con- 
tent of the so called theory of "inverse systems" introduced by Macaulay 
in |Macl916] . and developed in a context closer to our situation under the 
name of "Grobner duality " in [(^19701 [MMM19951 [EM2nn7j among others. 



The following is a multivariate analogue of Definition [27 

Definition 3.2. Let A := {Ai, A^} be a basis of A* as a K-vector 
space. Given any set E = {a;"^, . . . , cc""} of u monomials and given any 
polynomial h G -^^[aj], the generalized Vandermonde matrix Ve(A) and the 
generalized Wronskian matrix Wh,E{^) corresponding to E, A and h are 
the following u x D matrices: 



D 



Ve(A) = 



Ai(a;°i) • 




U 




Ai{x"^h) ■ 


■ Aoix'^'h) 


Ai(a;°") • 


■ Aoix'^-) 






Ai{x""h) ■ 


■ Aoix'^-h) 



We modify the definition of the matrix Os{Z) in (jl2p as follows: 
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Definition 3.3. Let S = {x'^^ , . . . ,x'^''} C K[x]t be of cardinality k 
T~iDi...Dn+iit), T* := Uj^^^-^Tj as in ^ and IZ as in (fiO|). Then 



05(A) 



Ks(A) 
^r-(A) 



e K 



DxD 



Note that by our assumption on T being a basis of A and A being a 
basis of A*, we have det (V7-(A)) 7^ 0. The following is the extension of 
Theorem 13. II to the multiple roots case. 

Theorem 3.4. Let {fi, ■ ■ ■ - fn+i) C K[x] and T := ^j>oTj specified in 
(0j satisfying our assumptions, and A be an arbitrary basis of A*. For 
any t € Z>o and for any S = {x'^^, . . . , x'^''} C of cardinality k = 

'HDi...D„+iit), we have 



^sifi, fn+l) 




det (05(A)) 
det {Vr{A)) ' 



Proof of Theorem \3.4\ The proof is similar to the proof of Theorem 3.2 in 
|DKS2006] . to which we refer for notations and details. Extra care must be 
taken however, as we are not anymore considering the polynomials /i, . . . , 
to have simple common roots. 

Using the exact same argument as in the proof of Theorem 3.2 in |DKS20d6] 
we can prove that 

± Sit) Asifi, . . . , fn+i) det (Vr(A)) = ± det(M') det (05(A)) , 

where 



(13) M' := 



is the submatrix of ([8]) obtained by removing the columns corresponding the 
monomials in T. In [DKS2006j we also showed that 

det(MO = ±£{t) ( H Ar^ 

^ j = t~Dn + l + l 

SO the claim is proved when £{t) 7^ 0. 
If £{t) = 0, we consider a perturbation "a la Canny" as in |Canl990] . i.e. 
we replace fi by fi^x := fi + \ x-' G i^(A)[a;], where A is a new param- 
eter, for 1 < i < n. It is easy to see that this perturbed system has no 
roots at infinity over the algebraic closure K{X) of K{X), since the lead- 
ing term in A of the resultant of its homogeneous components of degrees 
Di, . . . , Dn does not vanish, and hence the dimension of the quotient ring 
Ax := K(X)[x]/{fi^x, • • • ) fn,\) as a K{X)-vector space is also equal to D. 
It can also be shown (see [Canl990] ) that £x{t) 7^ 0, where £x{t) denotes 
the extraneous factor in Macaulay's formulation applied to the polynomials 
fi,\, 1 < i < n. Indeed, if Ef is the matrix whose determinant gives £{t) 
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with rows and columns ordered properly, it is easy to see that the perturbed 
matrix is equal to Et + \ I, where / is the identity matrix. 
Therefore, the statement holds for this perturbed family: 



(14) A5(/i,A,...,/„,A,/n+l) 




det (^^(rA)) 

det {VriTx)) 



for any basis Ta of Al (here, A7: a = A^,^(/i a, • • • , /n,A))- 



The subresultants appearing in p4j) are polynomials in A, that, when eval- 
uated in A = 0, satisfy: 



^s{fl,\, ■ ■ ■ , fn,\, /«+i)Ia=0 ~ ^sifl, ■ ■ ■ , fn, fn+l), ^Tj,X 



= Ar,, Vi. 

A=0 ^ 

So, in order to prove the claim, it is enough to show that there exists a basis 
of ^4^ which "specializes" to A when setting A = 0, i.e. to find a basis Aa 
of A*^ such that 

(15) det(05(AA)) 



det (05(A)) 



det (W(A)) ' 



det (Fr(AA)) 

and then to apply Identity to Aa and to specialize it at A = 0. 

We now construct the basis Aa: The monomial basis T = {a;°^, . . . ,a;"^} 
of A is also a monomial basis of Ax, since clearly linearly independent, and 

therefore it defines the dual bases {Vai, • • • > Vao} ^* {?/ai,A' • • • > 2/a_D,A} ^ 
A*^, satisfying for 1 < j,k < D, 

y^^{x"^) = 2/c»fc,A(^"^) = 1 if ^ = J arid otherwise . 



We write Aj = 'Ylik=i '^ikUaj, for 1 < i < D, where Cik G K, and then set 
Aa := {Ai,A, • • • , Ad,a}, with A^^a := J2k=i ^ikVck^^ l<i<D. Note that 

(16) K,x{x"^) = Ai{x"^) = c^j for l<i,j <D, 

(17) det (Fr(AA)) = det {VriA)) = det G K \ {0}, 

as the matrix {cij)i<ij<D is invertible. This implies that Aa is a basis of 

A* ~ ~ 

^A- 

We claim now that, for every a S N", there exist polynomials pa and Aj^a., 
I <j < D,m K[X] such that 

D 

(18) Pe,(A)x" = ^A,-„(A)a;°^ in ^A, and p„(0) 0. 

i=i 

For this, it suffices to express the monomial x" in terms of the basis T of 
Ax and take Pa{X) as a common denominator when lifting the expression to 
-fC[A][a;], satisfying the condition gcd (p^, Aj^cn 1 < J < ^) = 1- It is clear 
that Pct(O) 7^ because T is also a basis of A, and by assumption is not a 
common root of Pq. and Aj^a, ^ ^ 3 ^ D. 

Applying Aj^A to Identity (fT8]l and Aj to Identity (fT8|) specialized at A = 0, 
we then get by ()16p for 1 < i < D: 

D D 

Po.{\)K,\{x") = ^CijAj^oc{\) and p«(0)Ai(a;") = ^Cij^j>(0). 
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This imphes that the entries of the matrix Os{-^x) are the same JC-hnear 
combinations of the quotients "^^'"^^^ than the entries of the matrix 05(A) 

in terms of "^^^gy^- This and Identity (|17p imphes (jlSp . which proves the 
statement. □ 

As we mentioned before, for an arbitrary basis A of A* the expression in 
Theorem 13.41 may not provide a formula in terms of the roots of /i, . . . , /„. 
In order to obtain one, we recall here the notion of Grobner duality from 
IMMM1995] . 

For a = (ai, . . . , an) € N" define the differential operator 

On 



ai! dx^^---dx^" 

and consider the ring := {X^Qg^n aada : € K}. 

For 1 < i < n define the JT-linear map 



ai : K[[d]] ^ K[[d]]; aiid. 



d{ai,...,a,-l,...,an) if > 0, 

otherwise, 



and for /3 = . . . , G N" define = cif^ o • • • o a^. 

A K-vector space V C is closed if dimi^(y) is finite and for all (3 £ N" 

and D G F we have o"^(D) G V. Note that K[[c?]] and its closed sub- 
spaces have a natural i('[a3] -module structure given by x^'D{f) := Y){x^f) = 

Let ^ G K^. For a closed subspace V C define 

V€(F) := {/ G K[x] : D(/)(^) =0, VD G F} C K[x]. 

Let C K[x] be the maximal ideal defining ^. For an ideal J C define 

A^{J) := {B e K[[d]] : D(/)(0 = 0, V/ G J} C ^[[5]]. 

Then the following theorem gives the so called Grobner duality: 

Theorem 3.5 ( |Grl9701 [MMM1995] ). Fix ^ G K". The correspondences 
between closed subspaces V C ^i^[[i9]] and m^-primary ideals Q, V 1-^ V^{V) 
andQ^ A^{Q) are 1-1 and satisfy V = A^(V^(y)) andQ = (V^(A^(Q))). 
Moreover, 

dimi^(A^(Q)) = mult(Q) and mult(V^(y)) = dimi<-(y). 

We set Z := {^^, . . . , ^^j} for the set of all common roots of /i, ...,/„ in K^, 
with multiplicities, and C K[x] for the maximal ideal corresponding to 
for 1 < i < m. 

Example 3.6. f |EM2007l Exemple 7.37]) 

Let fi = 2x1X2 + 5x^, /2 = 2x\x2 + 5x| G C[x, y\. 

Then Z = {0, ^j^, ^3; ^4' ^5} where = (0,0) has multiplicity eleven and 
= (5|i")5|^) where ^ is a primitive 5-th root of unity, are all simple, 
1 < i < 5. 

Denote by Qo, Q^^, 1 < ^ < 5, the primary ideals corresponding to the roots, 
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then A^.((5^,) = (1) for 1 < i < 5 and if {61,62} is the canonical basis of 

7/2 



(494ei ~ 59ei+2e2)) (494e2 ~ 5(92ei+e2)) (3c^2ei+3e2 — C?5ei — d^e-})) ■ 

Using Grobner duahty, we are now able to give an expression for the 
subresultant in terms of the roots of /i, . . . , /„. For D G i^[[5]]and ^ € -ftT", 
we denote by D|^ the element of defined as D|^(/) = D(/)(^). In 
particular, under this notation, 1|^ = ev^. 

Corollary 3.7. Using our previous assumptions, let I = (/i, . . . , fn) and 

I = Qin---nQm 

be the primary decomposition of I, where Qi is a -primary ideal with 
di := mult(Qi). For 1 < i < m let Vi := A^^{Qi) c'K[[d]] be the corre- 
sponding closed subspace, and fix a basis {Dj^i, . . . , Dj^^^-} for Vi such that 
Dj 1 = 1. Then 



A :- {Di^il^^, . . . ■ ■ 

is a basis of A* over K. 



...,D 



m,dr, 



km} 



Note that the above choice for the dual basis A contains the evaluation 
maps for the roots of /, and using this A in Theorem 13.41 gives an expression 
for the subresultant in terms of the roots of /. 

Example 3.8. This is a very simple example containing an expression for 
a subresultant in terms of the roots. 

Let /i := X1X2, /2 := xf + (x2 - 1)^ - 1, /s := cq + ciXi + C2X2, with 
Co, ci, C2 € C. Then Z = {(0,0), (0,2)} where (0,0) has multiplicity 3 and 
(0,2) is simple. By computing explicitly, we check that T := {1, xi, 0:2, x|} 
is a basis of A = C[xi, X2]/{fi, f2) and that 

A := {l|(0,0), f^ei 1(0,0)) (f^e2 + 292ei ) | (0,0) , 1 1 (0,2) } 

is a basis of A* . We will use these bases to express the degree t = p = 2 
subresultant A^2(/i, /2, /s) with S = {xf} in terms of the roots of /i,/2. 
First, A^.2(/i, /2, /a) is equal to the following 6x6 determinant (since here 
the extraneous factor is 1): 



A,2(/i,/2,/3) = detM5 = det 












1 














-2 


1 





1 














1 





Co 


Cl 


C2 














Co 





Cl 


C2 











Co 





Cl 


C2 



C^ + 2C^C2. 
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On the other hand, Theorem 13.41 gives the following expression: 




det 

det 05(A) 


2 

Co Ci C2 Co + 2C2 

Co 2ci 
Co 2co + 4c2 


4(c[l + 2c2c2) 


det yr(A) 

det 


10 1 
10 
12 
4 


4 



using that Ti = is the degree 2 part of T, and 



Ar2(/i,/2)=det 











1 


det 


1 





1 







1 
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